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Abstract. Interferometry with Heisenberg limited phase resolution may play an 
important role in the next generation of atomic clocks, gravitational wave detectors, 
and in quantum information science. For experimental implementations the robustness 
of the phase resolution is crucial since any experimental realization will be subject to 
imperfections. In this article we study the robustness of phase reconstruction with 
two number states as input subject to fluctuations in the state preparation. We find 
that the phase resolution is insensitive to fluctuations in the total number of particles 
and robust against noise in the number difference at the input. The phase resolution 
depends on the uncertainty in the number difference in a universal way that has a 
clear physical interpretation: Fundamental noise due to the Heisenberg limit and 
noise due to state preparation imperfection contribute essentially independently to 
the total uncertainty in the phase. For number difference uncertainties less than one 
the first noise source is dominant and the phase resolution is essentially Heisenberg 
limited. For number difference uncertainties greater than one the noise due to state 
preparation imperfection is dominant and the phase resolution deteriorates linearly 
with the number difference uncertainty. 



1. Introduction 



The precision with which phase shifts can be determined in an interferometer is 
ultimately limited by shot noise. The shot noise originates from the discreteness of 
the interfering particles For interferometry with independent particles, each particle 
interferes only with itself. For N independent particles the precision is \/N times 
better than the single particle precision, a result which follows from application of the 
central limit theorem. This scaling of the phase precision with the number of particles, 
^0SQL ~ 1/V^5 is called the standard quantum limit. 

In some applications an increase in the number of interfering particles is difficult 
or impractical. In these cases it is worth trying to make better use of the available 

I For conciseness we refer to the interfering objects as particles with the understanding that they could 
be realized as actual particles or field quanta such as photons. 
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resources. Interferometry with phase resolution better than the standard quantum hmit 
is relevant for gravitational wave detectors [H [21 13|. The best atomic clocks employing 
Cs in E], trapped ions [HI [7], and earth-alkaline atoms such as ^''Sr in optical lattices 
[HI [9l [TOl [m [12] are also operating very close to the standard quantum limit. In the 
next generation of clocks the precision could be further increased by using interferometry 
with resolution better than the standard-quantum-limit. 

Phase resolution better than the standard quantum limit can be achieved with 
entangled many particle states. A variety of such input states have been proposed. For 
interferometry with photons, squeezed states of light have been suggested. Bollinger 
and Wineland have proposed spin squeezed states for ions and other massive particles 
|13j . Maximally entangled states have been proposed theoretically [lH [15] and their 
improved phase resolution has been demonstrated experimentally [T6]. Holland and 
Burnett have suggested to use Fock states in each of the input ports and they have 
shown that these states yield Heisenberg limited phase resolution [17]. States with 
improved phase resolution have been systematically studied in [18], [191 [201 [211 [22] ■ 

For linear interferometers the resolution achievable with N particles is 
fundamentally limited to a scaling A0h.l. ~ 1/A^, called the Heisenberg limit. 

So far, most studies of interferometry with resolution better than the standard 
quantum limit have concentrated on ideal realizations in which the initial state can be 
realized with perfect fidelity, no dephasing or decoherence happens during the evolution 
through the interferometer, and the final state is detected with ideal detectors. For 
the maximally entangled Schrodinger cat state, the effect of dephasing during the 
evolution through the interferometer has been studied by Huelga et al [23]. The effect 
of decorrelation on the phase resolution for interferometry with balanced Fock states 
has been studied by Kim et al [21] and recently Dorner et al |25j have constructed 
generalizations of maximally entangled states that are more robust against photon loss 
inside the interferometer. 

The goal of this article is to study the robustness of the dual Fock state against 
imperfections in the creation of the initial state. We find that this state is insensitive 
to fluctuations in the total number of particles fed into the interferometer. The phase 
resolution is robust against fluctuations in the number difference. If the fluctuations in 
the number difference are of order one, the phase resolution is still essentially Heisenberg 
limited. For increasing imbalance uncertainty, the phase resolution degrades linearly 
and becomes standard quantum limited as the imbalance uncertainty approaches the 
partition noise. These findings add a whole class of states to the arsenal of states 
that are candidates for phase resolution better than the standard quantum limit. This 
class of states can be thought of as the incoherent cousin of the conventional squeezed 
states. The states are incoherent mixtures of number difference states with a probability 
distribution that is narrower than the partition noise. 
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Figure 1. Schematic of a Macli-Zetinder interferometer. Particles prepared in state 
IV'in) pass through a beam splitter. The particles going through one arm acquire a 
phase shift (f) relative to the other arm. After recombination at the second beam 
splitter the particles coming out of ports one and two are detected. 

2. Model 

For the remainder of this article we consider the elementary Mach-Zehnder 
interferometer depicted in figure [H Particles prepared in state pass through a 
beam splitter. The component passing through one of the arms experiences a phase 
shift relative to the other arm. The two modes are recombined at a second beam 
splitter and detectors measure the numbers of particles coming out in ports one and 
two. 

We describe the particles passing through the interferometer as effective two level 
systems with the two states corresponding to the two arms of the interferometer. 
Mathematically we describe these two level systems using the angular momentum 
formalism. This formalism as it applies to interferometry has been described in detail 
in several articles, see for instance [21], [211 126]. In brief, each particle is identified with 
a fictitious spin-1/2 system with spin up and down corresponding to the particle being 
in port one and two, respectively. A collection of particles is then described by 
combining the N spins into a total angular momentum. 

The 2-component Jz of the total angular momentum describes the number 
imbalance between the ports of the interferometer and the square of the length of the 
angular momentum vector is related to the total number of interfering particles. We 
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choose the and i 

Jz\J,m) 



eigenstates \J,m) with 
> = m\J, m) 



and 



J^\J,m) 



J(J + 1)|J,M) 



and J = N/2 as basis states. The perfectly balanced Fock state suggested by Holland 
and Burnett is the state \J,0) The interferometer transformation that relates the 
output state to the input state is 



From here on we specialize to the case of an interferometer with relative phase shift 
= 0. For the balanced input states that we consider below this choice of working 
point corresponds to the null fringe were equal numbers of particles are detected in 
both output ports of the interferometer. 

2.1. State preparation 

In an ideal experiment it would be possible to deterministically create the state \J,0) 
with fixed J. However, for systems with large numbers of particles, where the precision 
gain would be most dramatic, this perfect situation cannot be achieved in general. 
In practice, states \J,m) are created stochastically with a probability distribution Vj^m- 
Typically it is impossible to determine which J and m were realized due to finite detector 
sensitivity. This means that one has to reconstruct the phase based on incomplete 
knowledge of the initial state. 

It has been shown by Uys and Meystre [2lj that the phase resolution for states | J, m) 
depends very little on m as long as |m| -C J. [Jj] We can therefore assume, without loss 
of generality, that the mean of the distribution Vj^m is centered at m = in accordance 
with the original proposal of Holland and Burnett [I7j. For the working point = 
this is also the most interesting case from an experimental perspective, where it is 
advantageous to detect approximately equal numbers of particles in both detectors due 
to finite detector efficiency. Modeling the distribution Vj^m as a Gaussian distribution 
in J and m we can characterize the state preparation by three parameters: the average 
total angular momentum J and the standard deviations AJ = "Pj'iJ' ~ JY 

and Am = a/^^ Vm^"^ of the distributions of the total angular momentum J, 
Vj = J2m'^J,m, and of the number imbalance m, Vm = J2j'^J,m- From now on we 
refer to Am as the state preparation imperfection or just imperfection for brevity. The 
initial state is then a statistical mixture. 



§ For simplicity we restrict ourselves to even total numbers of particles. The case of odd numbers of 
particles is qualitatively identical. 

II In fact, these authors have shown that the phase resolution is nearly Heisenberg limited almost all 
the way to |rn| J. 



|V'out)=e*^^-|V'in). 




(1) 
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2.2. Phase reconstruction: Bayesian inference and likelihood function 

It is well known that reconstructing the phase shift from a measurement record is a non- 
trivial task if the input state of the interferometer is the balanced Fock state. For that 
state, the mean number imbalance ( J^) is zero for all phase shifts and hence cannot be 
used as a signal. The variance (AJ^)^ varies sinusoidally with the phase, but the signal- 
to-noise ratio of that observable is of order one. Thus a large number of measurements 
are necessary in order to infer the phase shift. It can be shown that the same holds 
true for all higher moments of Jz as well. In order to achieve that state's potential 
phase resolution, one has to use a Bayesian reconstruction scheme instead [17]. This 
scheme has the additional advantage that it yields the full probability distribution of 
the phase conditioned on a measurement record instead of e.g. the first few moments. 
Furthermore, it immediately yields a prescription for how to find the phase in an actual 
experiment. 

The basic idea of Bayesian inference is to update the probability distribution of the 
phase according to 

with M a normalization constant. Note that we use the symbol M to denote different 
normalization constants throughout this article whose precise value is irrelevant. The 
prior phase distribution P(0) contains all the knowledge available about the phase before 
the first measurement. To provide an unbiased analysis of the phase measurement we 
assume that nothing is known about the phase initially except that it lies in the interval 
[— 7r/2, 7r/2]. This latter restriction removes the irrelevant phase ambiguity where 
phases that differ from each other by an integer multiple of vr cannot be distinguished. 
The prior phase distribution is then 



TT 



-1 



< 7r/2 



P(m|0) is the conditional probability for a number difference m to occur given that 
the phase shift was (j). The denominator normalizes the probability distribution P{(f)\m). 

Equation ([2]) can be iterated: After a measurement of the number difference mi 
the knowledge about the phase is P(0|mi). The phase distribution after a second 
measurement is then obtained from equation ([2]) by using P((/)|mi) as the prior. For a 
measurement record {mi,m2, ■ ■ ■} one finds the phase distribution 

p/,| ^ P(mi|0)P(m2|0) ■ ■ ■■ , . 

P(0|mi,m2,...) = , (4) 

The sharpness of this distribution after a certain number of measurements is a figure of 
merit of the phase resolution obtainable with a given input state. 

% This choice of the prior phase distribution makes this phase inference method essentiahy a maximum 
Hkehhood estimation. 
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In the limit of a large number of measurements n the various possible measurement 
outcomes m occur nP{m\6) times, where 6 is the actual phase shift. The phase 
distribution then becomes 

J 

P{(j)\e)=^-' Y[ P(m|0)"^('"l^) (5) 

m=—J 

= Ar-'F{<p\er, (6) 

with 

J 

mo) = n ^("^i^)''^'""^- (7) 

m=—J 

The function F(0|6') is the asymptotic likelihood function of the phase. It fully 
characterizes the expected phase resolution in the limit of large n. 



3. Results 



3.1. Conditional probabilities Pp{m\(j)) 

From equation ([2]) it is clear that the conditional probabilities P(m|0) are of central 
importance for the phase resolution. The sharper the distribution P(m|0) is as a 
function of 0, the more one learns about the phase if the number difference m is detected. 
However it is not sufficient for a state to have sharply peaked P(m|0) distributions 
in order to be useful for interferometry. Measurement outcomes m with such sharp 
distributions also must be likely to occur. This latter requirement makes the phase 
resolution typically dependent on the working point of the interferometer, i.e. the 
precision depends on the phase that one wishes to measure. 

The conditional probabilities P(m|0) are calculated by sending the initial state p of 
equation ([1]) through the interferometer and projecting the final state on the subspace 
with a given m, 

Pp{m\<f)) = Tr (e^^-^-pe-^-^-^-n^) , (8) 

where 11^ is the projector on the m subspace, 

ftm = ^ \J,m){J,m\. 
J 

The conditional probabilities for a well defined initial J and m, P\j,mo){J,mo\{''T^\4') , behave 
as 

P\J,mo){J,mo\M(P) = Jlo-miJ^P) (9) 

for phase shifts small enough so that the number differences m near the poles of the 
Bloch sphere are unlikely. In this equation Jm^-m is the Bessel function of the first kind 
of order mo — m. 

In Figure [2] we show the conditional probabilities for various initial states for an 
average of 20 particles, i.e. J = 10. As a consequence of the Bessel function like 
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Figure 2. Conditional probabilities Pp{m\<j)) for initial states p for different levels of 
uncertainties A J and Am for J = 10. In the first row there is no uncertainty in the 
total particle number, A J — 0, and the number difference uncertainties are (a)Am = 
(perfect state preparation), (b) Am = 1, and (c) Am — 3. With increasing state 
preparation imperfection Am the central peak of the distribution becomes broader 
leading to reduced phase resolution. Panels (d-f) show the probability distributions 
for the same state preparation imperfections as the panels in the first row but with 
AJ = 3. These panels demonstrate that additional uncertainty in J does not affect 
the central peak of the probability distributions. Rather such noise affects the wings 
of the distributions for phase shifts (p ^ 1/A J. 

behavior, the conditional probabihty P(O|0) for perfect state preparation in figure [2]^a) 
is sharply peaked at = with a width of order ~ 1/J. This sharp peak gives rise 
to Heisenberg limited phase resolution for the perfectly balanced Fock state. For larger 
phase shifts the conditional probabilities oscillate with frequency J and with an envelope 
falling off as 2/(7rJ0). 

Figures [2](b,c) show that the width of the peak of the conditional probabilities 
increases with increasing state preparation imperfection Am. Qualitatively, it is clear 
that this increased width leads to a decreasing phase resolution. The width of the 
central peak can be estimated by the construction illustrated in figure [3]^a). For phase 
shifts 101 <^ 7r/2 the maxima of the distribution P{m\(j)) follow the red lines given by 
m = ±J(j). Mathematically, this follows from properties of the Bessel-functions which 
describe the conditional probabilities well for |0| ^ 7r/2. Qualitatively, the scaling of m 
with can be understood by noting that the number difference has to increase from zero 
to \m\ = J for a phase shift of ±7r/2. These lines translate the Gaussian distribution 
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Figure 3. Figure (a) illustrates the construction of the estimate of the width of 
Pp{0\(f>) for J = 20, AJ = and Am = 3. The red lines indicate curves m = ztJcj). 
The two arrows indicate the widths of Gaussian distributions along m and (j). Figure 
(b) shows the agreement between the central portion of Pp{0\(f>) and a Gaussian with 
width Acf) = Am/J normalized to the same height (blue solid line) for parameters 
identical to (a). For phase shifts larger than 1/J good agreement is found between 
Pp{0\(j)) and {nJ\(j)\)~^ (purple dashed line) which follows from the asymptotic behavior 
of the Bessel functions. 



PpH<P) - zm- (11) 



P(m|0) of m with width Am into a Gaussian distribution of the phase with width 
Am 

A0 = (10) 

This Gaussian is compared with the exact probabihty distribution P(O|0) in figure [3]^b). 
The central portion of that probabihty distribution is weU described by the Gaussian. 
The wings of the conditional probability distributions can be found from the asymptotic 
behavior of the Bessel functions. Averaging several Bessel functions leads to the tails 
falling off as 

1 

nJcj) 

This is also illustrated in figure [3]^b). Note that in a phase reconstruction experiment 
the wings get more and more suppressed with increasing number of measurements and 
eventually the phase distribution is entirely determined by the approximately Gaussian 
central peak. 

The figures [2](d-f) show the same conditional probabilities as in the panels (a-c) but 
with AJ = 3. As a consequence of the addition of Bessel functions that asymptotically 
oscillate with different frequencies the oscillations in the tails wash out for phase shifts 
101 > 1/AJ in figure [2](d). Near the poles |m| ~ J of the Bloch sphere the distributions 
look different compared to case with a sharp total number of particles, A J = 0, because 
the total angular momentum now has a spread. Most importantly however, the central 
Gaussian peak is unaffected by the noise in J. This is true for noise up to AJ ~ J. 
As a consequence, the phase reconstruction is largely insensitive to this type of noise. 
Physically, this is because the measurement is based on detection of number differences 
in which the total number of particles cancels out. Only when the phase shift is so large 
that almost all particles exit the interferometer in one port does the total number of 
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particles play a role. We have numerically confirmed this independence of the phase 
resolution of AJ. That allows us to consider the limit AJ from now on without 
loss of generality. 

3.2. Simulation of phase reconstruction: Finite numbers of measurements 

To illustrate one of the key new features of phase reconstruction with imperfect state 
preparation compared to the ideal case we simulate an actual phase reconstruction 
experiment. As in the numerical examples in the previous section we consider a total 
particle number of iV = 20 and the state preparation imperfection is Am = 1. 

To simulate the phase reconstruction we draw a sequence of measurement results 
m at random with probabilities given by Pp{m\(j) = 0). Using the known conditional 
probabilities Pp{m\(j)) we then successively update the phase distribution by means of 
Bayes' formula. The evolution of the phase distribution with successive measurements 
is not deterministic due to the stochastic process with which the measurement results 
are determined. 

Two representative examples of such phase reconstruction simulations are shown 
in figure (jl]). The first example is reminiscent of phase reconstruction with perfect 
state preparation. The phase distribution is centered at the correct phase shift = 
and the width of the distribution shrinks monotonically with increasing number of 
measurements. 

The second example illustrates a complication that can happen with imperfect 
state preparation. In the fourth measurement of that sequence a number difference of 
m = —2 is detected. For the state preparation imperfection Am = 1 this measurement 
outcome is relatively unlikely to occur for zero phase shift and the Bayesian updating 
wrongly attributes the measurement result to a non-zero phase shift. The phase 
distribution becomes wider and is peaked at non-zero phase shift. Thus this unfavorable 
measurement outcome reduces the phase resolution and, even worse, it leads to an 
inaccurate phase inference. Further measurements correct this measurement and in the 
limit of large numbers of measurements the phase inference is accurate. 

General quantitative statements on how many measurements are necessary in order 
to rule out such inaccuracies for a certain level of state preparation imperfection are 
difficult to make. In practice it is probably easiest to simply run numerical simulations 
of the phase reconstruction in order to empirically decide how many measurements 
should be made. 

3.3. Asymptotic phase resolution 

Asymptotically, the phase resolution A^ improves as ~ Xj^fn^ where n is the number 
of measurements. A figure of merit that is independent of the number of measurements 
and characterizes the suitability of an input state for phase measurements is therefore 




(12) 



■oo 
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0.5 0.5 

Figure 4. Two different runs of numerical phase reconstruction simulations. The 
figures show the phase distributions P{4)\mi, m2, . . . , as a function of the length 
of the measurement record, n. Panel (a) shows a reconstruction run that is similar 
to what one finds for interferometry with perfect state preparation. In panel (b) an 
unfavorable measurement outcome of m = — 2 in the fourth measurement leads to 
a wrong inference of a non-zero phase shift. This inaccuracy is remedied by further 
measurements. 



where A0(n) is the standard deviation of the phase distribution after n measurements. 
We can calculate A0oo through the variance of the phase distribution P{(j)\6) given in 
equation with an n that is large enough so that the scaling equation (fT2l) holds. For 
the results that follow we have used n = 5. This asymptotic phase distribution is devoid 
of the statistical fluctuations encountered in the previous section for small numbers of 
measurements. 

We have numerically calculated the asymptotic phase resolution A(f)oD as a function 
of Am for different values of J. The results are shown in figure [5l In that figure we 
have rescaled the phase resolution by J~^. Remarkably, the resolutions for different 
J rescaled this way all collapse onto a single universal curve. For nearly perfect state 
preparation the resolution saturates at the Heisenberg limit, 

A0ooJ ^ 1, for Am < 1. (13) 

For Am > 1, A0oo increases linearly with Am. 

The universal behavior of JA(f)^ can be understood by treating the fundamental 
1/ J noise due to the Heisenberg limit and the noise due to the state preparation 
imperfection as two independent noise sources. The total noise is then obtained by 
adding these noise sources in quadrature. According to the construction explained in 
figure [3] the noise due to the state preparation imperfection is given by a Am/ J, where a 
is a constant of order one. a describes by how much F{(j)\6 = 0) is broadened compared 
to the Gaussian in figure Mjo) due to the admixture of results m 7^ 0. Adding then the 
Heisenberg limit noise and the state preparation noise in quadrature we find 

A0OO = V-J'^ + (aAm/Jf. (14) 

This scaling of the phase resolution with J and state preparation imperfection Am is 
the central result of this article. We find excellent agreement between equation f|T^ and 
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Figure 5. Numerically determined asymptotic phase resolution A(/)oo as a function of 
the state preparation imperfection Am. The figure shows the resolution rescaled by 
for J — 10 (blue discs), J = 20 (purple squares), J = 40 (yellow diamonds), and 
J = 80 (green triangles). The resolutions rescaled this way fall on a universal curve 
in agreement with formula (fT4|) . The inset shows the data points for Am = 0.8 at a 
higher resolution. The resolution described by equation (|14|) with a = 2 is shown by 
the purple line. 



the numerically determined phase resolutions for a = 2, as illustrated by the purple line 
in figure [H 

3.4- Phase reconstruction with unknown Am 

In a typical experiment one will not precisely know the initial number difference 
uncertainty Am. Especially in light of the numerical simulations for finite numbers 
of measurements presented in section 13.21 it is important to understand the dependence 
of the phase inference on the choice of an estimate for Am. In particular one needs 
to know under what circumstances the phase distribution is asymptotically accurate. 
Or, phrased differently, how should one estimate Am so that unfavorable measurement 
results like in figure |4](b) are eventually corrected by a large number of measurements. 

This problem can be studied by looking at the likelihood function F{(f)\6) of equation 
(I7j) as a function of the estimate Amest for the number difference uncertainty. Basically 
one pretends that the number difference fluctuations are Amgst and calculates the 
corresponding conditional probabilities -Pp(Amost)('^l0) using the initial state p{Amest) 
with imperfection Amgst- But in an actual experiment each measurement result occurs 
with probability Pp(Am){fn\4>) where the true fluctuations Am may differ from our 
estimate Amest- In the limit of a large number of measurements the phase distribution 
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Figure 6. Likelihood function F{(j)\d = 0; Amcst, Am) of the phase as a function of the 
estimate Amest for the number difference uncertainty for J = 20. The actual number 
difference uncertainty Am = 2 is marked by the black line. For too optimistic estimates 
Amcst < Am the likelihood function becomes bimodal and the asymptotic phase 
distribution is inaccurate. For conservative estimates Amcst > Am the asymptotic 
phase distribution is accurate but broader than for Amcst = Am. 

is then governed by the hkehhood function 

F{(P\9- Am,,,, Am) = A^-^ J] Pp(A„^..o(^l0)'''''^'"'^""^• (15) 

m 

That function is shown in figure ([6]) for J = 20, Am = 2, and ^ = 0. The figure 
illustrates two typical properties of F{(f)\9 = 0; Amest., Am). If the estimate of the 
fluctuations is too optimistic, i.e. if Amgst < Am, the hkelihood is peaked at a non-zero 
phase shift and it has a minimum at zero phase shift. For such estimates the phase 
distribution is asymptotically inaccurate. This is because for a too optimistic estimate 
of Am one observes non-zero number differences at the output too often to be brought 
about by a non-zero phase shift. One thus wrongly concludes that the phase shift was 
non-zero. 

For conservative estimates Amest > Am on the other hand the likelihood function 
of the phase is peaked at zero but it is broadened compared to the likelihood function 
for perfect knowledge of the number fluctuations, Amest = Am. As expected, the best 
phase resolution is obtained for Amcst = Am. 



Robustness of Heisenberg-Umited interferometry with balanced Fock states 



13 



4. Conclusion 

Interferometry with balanced Fock states is robust against imperfections in the state 
preparation. The phase resolution is insensitive to noise in the total number of particles. 
For number difference uncertainties smaller than one the phase resolution is essentially 
Heisenberg limited. For larger state preparation imperfections the phase resolution 
deteriorates linearly with the number difference uncertainty. 

In future work we plan to study the robustness of interferometry with balanced 
Fock states against dephasing and decoherence inside the interferometer as well as the 
effect of non-ideal detectors. 
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